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I. introduction 


Background « The genealogy of the work reported in this paper stems 
From a mixing zone research project originated in 1968 and supported 
financially by the Department of Natural Resources of the State of 
Wisconsin. The central objective was the study of the character- 
istics of dispersion and diffusion as the mechanisms by v/hich pol- 
lutants are transported in natural river courses, with the view of 
providing additional data for the establishment of water quality 
guidelines and effluent outfall design protocols. 

The first year involved a comprehensive study of the literature 
and some laboratory model observations from v;hich there developed a 
compendium of reference materials, some qualitative conclusions 
regarding the effects of certain outfall, effluent and river char- 
acteristics on mixing rates, and a functional hypothesis for the 
extent of the mixing zone in terms of these characteristics. The 
second year included the first field surveys (during the summer of 
1969) together with the beginnings of laboratory model studies of 
one of the field sites. At this time also the mixing zone project 
became part of a larger interdisciplinary Remote Sensing Program 
administered by the University’s Institute for Environmental Studies 
and financed by the Office of University Affairs of the National 
Aeronautics and Space Administration (NASA) . This evolvement was 
natural, because the summer field surveys of 1969 brought out the 
inadequacy of ground activities alone for covering large geographic 
areas, namely that observed data are largely either static or time 
averaged- Over the last three years additional ground and aerial 
field surveys and laboratory modelling have been unaortaken. in 
addition, mathematical modelling and quantitative photogx<xnhic 
delineation of the mixing zone have been initiated. 

Project Goals . Throughout this investigation, work has been divi^g^ 
into four basic categories which are directed at the basic goal ot 
developing relationships which will permit the estimation of the 
nature and extent of the mixing zone as a function of those variables 
which characterize the outfall structure, the effluent, and the 
river, as well as climatological conditions. These relationships 
may be used (a) in the establishment of definite and rational water 
quality guidelines, (b) in the development of a sampling and surveil- 
lance program by governmental and private agencies, and (c) in the 
development of a protocol for the design and location of waste 
effluent outfalls. 

The four basic categories of effort are; 

1. The development of mathematical models. 

2, Laboratory studies of physical models. 
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3. Fi.eld surveys involving ground and aerial sensing. 

4. Correlation between aerial photographic imagery and mixing 
zone characteristics. 


In category (1) simple jet- and source-type, two- and three- 
dimensional mathematical models have been constructed and compared 
with field observations. Detailed mathematical modelling, including 
more realistic approximations of the flow field and river and 
outfall geometry, is being undertaken. The work in Category (2) 
has led to observations and measurements of temperature patterns 
in laboratory models of the Vvfeston power plant (see site descrip- 

« — L^nS/ diiu wl. ^11 ou^ii TfiS WOjTjC 

in Category (3) includes field measarements of effluent discharge 
mixing patterns at seven sites throughout the State. These surveys, 
which have been carried out over the past four years, will be con- 
tinued with intensive ground^-aeriai measurements at one site. In 
Category (4) , a method has been developed vrhich relates aerial 
photographic film response (film density) to surficial suspended 
solids concentrations in the paper mill discharge at the Kimberly- 
Clark plant and the sewage treatment plant (STP) discharge at 
Neenah-Menasha (see site descriptions) . The work in this category 
will be intensified and move in somevjhat different directions as a 
result of the work done tlras far. 


Scope of this Report . The general objective of this report is to 
give an overview of the pirojecrt goals and activities to date, 
together with a conceptual dJ-scussion of the mixing zone definition 
and measurements associated r,herewith. Accordingly it will be 
necessary to shov; the powex'fully good consequences of coordinating 
ground and aerial techniques and methods, in contrast to their 
inadequacies v/hen separated* Xn addition, sufficient results will 
be presented to serve as indicators of the type of v7ork accomplished 
and yet to evolve as the work proceeds. Three of the seven field 
sites will be used in order to cover a fairly v/ide spectrum of waste 
effluent situations indigenous to the industries and municipalities 
of the State of Wisconsin (see site descriptions) . 


II. DEFINI TION OF THE MIXING ZONE 


T<<rhenever two or more fluid substances are brought together, a 
mixing zone or region is created vrithin the regime of flow. The 
characteristics of the mixing zone are related to the convective 
transport mechanisms of diffusion and dispersion. The transport 
associated with diffusion can be due either to molecular or turbulent 
actions, in most real life situations the "time" is so short that 
the action of molecular diffusion can be neglected. Therefore, in 
the mixing zone, diffusion is only the result of turbulent action. 
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If the fluids which come together are composed of identical 
substances, the characteristics of the mixing zone can be measured 
only by the momentum and energy distribution coefficients and the 
turbulent energy spectra. Manifestations of these are the velocity 
and energy gradients throughout the regime. Tvro homely examples 
of this situation are the confluence of rivers and the impingement 
of the air from a ventilating duct outlet into the air of a room, 
where the compositions of the fluids, in both examples, are the same. 

If the fluids which come together are not composed of identical 
substances, then the characteristics of the mixing zone can be more 
completely defined by the additional measurement of the distribution 
of the concentration of one of the substances within the other. 

The most common examples of this situation involve the introduction 
of effluent wastes into lakes, rivers and streams. The effluent 
wastes from domestic sewage treatment plants, industrial plant 
operations, and from thermal power generating plants are the most 
common in relation to their impacts on the earth's environment. 

Since the mixing zones of the second case can be measured more 
completely utilizing changes in concentration, it is practical to 
convert Case 1 into Case 2 through the introduction of tracer 
elements, which can be either chemical or particulate in nature. 

Thus it is possible to use fluorescent dyes and nuclear isotopes in 
such minute quantities that they do not affect the diffusion and 
dispersion transport characteristics, but their concentrations can 
be measured easily with adequate accuracy. 

The physical mixing zone region is defined as that portion of 
the flow regime in which the tracer (or effluent waste) is diluted 
to the concentration obtainable by total or complete mixing at the 
point of discharge. Such a definition presupposes that the flow 
regime and the tracer regime have the initial properties of real 
fluids in motion (at least mass density and absolute viscosity) and 
are subjected to the constraints of real physical boundaries. 

The dilution processes are related to the convective transport 
mechanisms of diffusion and dispersion, as mentioned earlier. This 
mixing activity is accomplished in two stages, each having its own 
logical part of the mixing zone region. 

Stage one is always the first part of the zone. Here the 
exchange of momenta is on a macro-scale, since the initial conditions 
of velocity, linear momentum, turbulence levels, and turbulent 
energy spectra of the flow regime and effluent wastes are not the 
same (if they were, stage one would not exist) . Thus, the mixing 
and dilution processes are largely dispersive in nature, because the 
scattering is clearly visible and is manifest by macroscopic eddies 
and swirls clearly different in character than those of the aunbient 
regime. 
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Staqe two begins at the completion of the macro-scale activity, 
therefore the velocity, momentum, and turbulence characteristics of 
the zone and the ambient regime are about the same. They can't be 
exactly the same, however, since the two liquids ordinarily have 
different suspended solids, temperatures, and kinematic viscosities, 
each of which has an effect on the turbulence level for a given 
energy gradient. It naturally follows that the dilution process in 
stage two is on a relatively small scale and is the result largely 
of turbulent diffusion* 

Stage two of the zone ends v;hen all of the characteristics of 
the zone and the ambient are the same, or as stated previously, when 
the tracer is diluted to the concentration obtainable by complete 
mixing. It can also be said that stage two cannot be seen by the 
human eye, unless of course the effluent waste has visible charac- 
teristics different from those of the ambient regime - 


III. MEASURING THE MIXING ZONE 


The fluid dynamic, thermal, and water quality characteristics 
of the mixing zone can be quantified using both ground and aerial 
sensing methods in complimentary ways. By so doing it is important 
to realize the separate and coordinate roles which each plays. 

The ground measurements, although slow and tedious, define the 
vertical variations and provide the data for correlating aerial 
photographic imagery and effluent wastes concentrations- The aerial 
sensing methods relate to the collapse of the time scale* In a 
matter of a few minutes the surficial layer several miles long and 
wide can be photographed. In a sense the aerial measurements "stop 
the clock" and thus show the extensiveness of the zone, its tran- 
sient nature, as well as provide relative film densities. This is 
an excellent example of the benefits that can be derived when several 
disciplines are brought together in such a way as to maximize the 
results. 

In the presentation of results more detailed descriptions of 
the ground and aerial sensing techniques and methods of analysis 
are given. 


IV. DESCRIPTION OF FIELD SITES 

The Weston Power Plant . The plant is a 135 megawatt steam electric, 
fossil-fueled, unit located on a fairly straight reach of the 
Wisconsin River a few miles south of Wausau near Rothschild. The 
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plant is 40 percent efficient and takes in condenser cooling water 
from the river through an intake of traditional design, except that 
the water passes over drum screens to remove the suspended solids 
introduced by paper mill operations upstream near Wausau. The 
screenings are continuously flushed off the drxam and reintroduced 
into the river in the locale of the effluent outfall, as shown in 
figure 1. 

The v/arm water effluent is discharged through a reinforced con- 
crete structure, whose opening is 6 ft by 8 ft and submerged about 
3 ft, into a small bay. The bay is connected to the river by a 
canal whose centerline is pointing about 45 degrees doT^stream from 
a normal to the river flow pattern. 

The river width is about 500 ft, and under average conditions, 
the depth is about 5 ft, thereby producing a width to depth ratio 
of 100. The normal summer discharge of the river is 2000 cfs, with 
a velocity of about 1.0 fps. The discharge and velocity ratios 
(river to effluent) vrere 10 and 0.25, respectively. 

The Waukesha Sewage Treatment Plant . The municipal sewage from the 
City or Waukesha, Wisconsin, is handled in a plant having secondary 
treatment and located on the Illinois Fox River just south of the 
city. The effluent is discharged into the river on the left bank 
through a 30 -in. diameter culvert pipe aligned about 30 degrees to 
the axis of the river. In this region the river is very winding, 
and the outfall is located at the end of a large double bend, as is 
shown in figure 2. 

The river width is 80 ft on the average with a depth of 1.0 ft, 
so that the v/idth to depth ratio is 80. The normal summer 

river discharge is about 40 cfs, with a velocity of 0.7 fps. The 
discharge and velocity ratios (river to effluent) were 2.6 and 0.25, 
respectively . 

The Kimberly-Clark Paper Mill . The mill is located at the Cedars 
Dam on the Lov/er Fox River near Kimberly, Wisconsin. The plant 
manufactures special types of paper for color printing (photographic 
magazines, etc.) and discharges about 50,000 lbs per day of effluent 
wastes contining dyes, kaolins, titanium dioxide, pulp, etc. 

The mill has a total of 25 outfalls, four of which are located 
on the upstream side of Cedars Dam, and the remainder on the down- 
stream side. The discharges from the downstream outfalls combine 
into one plume along the south bank of the river, as shoxim in 
figure 3* About 65 percent of this discharge comes from the main 
mill outfall which is a U-shaped channel having an area of about 
3.2 sq ft. 

The river downstream from Cedars Dam bends to the right and 
widens from 600 ft to 1000 ft beyond the exit channel from the dam 



6 



This page is reproduced at the 
back of the report by a dilFerent 
reproduction method to provide 
belter detaiL 


Figure 1. 


Aerial Photo of Weston Power Plant Site. 
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Figure 2. Aerial Photo of Waukesha Sewage Treatment Plant Site. 
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Figure 3 . 


Aerial Photo of Kimberly-Clark Paper Mill Site 


9 


locks. About 0.75 miles downstream from the Cedars Dam is Little 
Chute Dam, which fact produces a mixing zone that enters the upstream 
or delta region of an impoundment. The river depth varies from 
5 to 10 ft. 

The normal summer river discharge is 1700 cfs with a velocity 
which is somewhat greater than 1,0 fps at the dam and less than 0.3 
fps upstream from the Little Chute Dam. The discharge and velocity 
ratios (river to effluent) were 120 and 0.20, respectively. 


V. PHOTOGRAPHIC AND THERMAL IMAGERY 


Background. Ground measurements at each of the sites were initially 
Supported by 35 mm color and color infrared photography obtained 
during exploratory aerial measurements taken in developing an inven- 
tory of pollution sites throughout the State. In most cases the 
photographs depicted in a rather vivid but qualitative way the 
presence, extent, and character of the mixing zone. 

Microdensitometric analysis of the photography of the Neenah- 
Menasha STP outfall suggested a very probable correlation between 
image densities and the level of concentration of suspended solids 
in the surficial layer (top 0.5 ft) of the mixing zone. It was 
decided to pursue the validity of the correlation by accelerated 
and more intensive measurements at the Kimberly-Clark Paper Mill 
site. 


In addition, aerial thermal scanning also has been utilized at 
several sites to extend and complement ground measurements of 
temperature patterns. The application and results of this method 
for the Kimberly-Clark Mill site are presented herein. 

Photographic Image Density I'lodel Development . Laboratory reflectance 
analyses performed by Scherz and Klooster(4) were used to demonstrate 
that water sample reflectance is a function of turbidity and the 
related suspended solids found in samples taken from the Kimberly- 
Clark mixing zone and ambient regions. From these data it can be 
shown that for a given wavelength and for solids concentrations 
less than about 200 mg/1, the relationship takes the form, 

R = K^S + K 2 (1) 

where R is the water sample reflectance, S is the concentration of 
suspended solids in mg/1, and and K 2 are parameters. 

The film exposure can be related to scene reflectance by the 
expression. 


E « K^R + 


(2) 
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where E is the exposure of the film, R is the scene reflectance, 
and and are also parameters* 

The image optical density resulting from relative film expo- 
sure can be determined on the basis of a film-step v/edge calibration. 
The resulting functional expression, as shown by Dana (2), can take 
the form, 

D = (3) 

where D is the image optical density? B^, B 2 , and B 3 are film 

type, handling, and processing parameters? and ^ is equal to rela- 
tive liOg E» 

Equations (1) through (3) permit the development of an expres- 
sion (model) which permits the estimation of suspended solids, S, 
from the measurable image density, D, so that, 

S = + 3 (4) 

where a and ^ are model parameters and ^(D) results from the solution 
of equation (3) . 

Application of Image Density Model . The model manifested by equa- 

T4) was used fn the study of the Kimberly-Clark Paper Mill 
mixing zone. The image optical densities were extracted from each 
of the three layers (red, green, and blue) of color infrared film 
type 8443. Concurrent with the aerial photography, suspended solids 
water samples were taken in the surficial layer at 82 locations. A 
photogramme trie solution for the image position associated with each 
sample location provided a viable means for obtaining simultaneous 
image densities and solids concentration. 

The image densities were digitally recorded using a scanning 
microdensitometer system, which permitted discrete density observa- 
tions over film areas depicting about 1,33 ft square areas on the 
ground. Experience showed that the image densities could be grouped 
into 36 density levels to give satisfactorily refined density 
observation (about 0.08 density units) . 

Equation (3) was solved for each observed density for corre- 
sponding values of Z, The B parameters were estimated from the D 
versus logE curves resulting from the film calibration v;ork. The 
suspended solids data resulting from the samples taken at the 82 
locations were then used with these ^ values to determine the a and 
3 parameters. In both instances above, the curve fitting was done 
by least squares methods. 

Figure 4 is a plot of the observed suspended solids and film 
density data. Also shown are the model predictions using a solution 
of equation (4) . it is believed that the scatter of the data is due 
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Fig. 4. Observed suspended solids versus relative image density 
for the Kimberly-Clark Paper Mill mixing zone. 
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largely to the experimental errors in obtaining and analyzing the 
water samples. 

The model (equation (4) or figure 4) can be used with confi- 
dence for this particular mixing zone, provided the a and 3 para- 
meters are determined for each roll of film. The applicability of 
this method to other types of suspended solids and mixing zones is 
a continuing objective of this program. 

Application of Thermal Scanning . Infrared thermal imagery obtained 
from thermal scanners caiTbe converted to surface water temperatures 
using the relation 

%r = A + BT + CT^ (5) 

where = scanner output (proportional to emitted radiation) , 

T = water surface temperature and A, B, C are parameters. Equation 
5 is a close approximation to the integral of Planck's lav/ for 
radiation over the 8—14 micron wave length range with corrections 
for the filter/ gain and offset of the scanner. To utilize equation 
5 and the 70 mm film output of the thermal scanner to determine 
water surface temperatures, the scanner film is first rectilinear- 
ized to correct for scanner distortions and enhanced to emphasize 
water temperatures and their variations relative to the land. Film 
density variations throughout the mixing zone are then determined 
using a microdensitometer and converted to a coded 256 level array 
(covering the maximum to minimum density variations) and stored on 
magnetic tape. Utilizing the known film density and water tempera- 
ture at several positions throughout the mixing zone, the parameters 
A, B, and C in equation 5 can be evaluated, Subsecjuently, water 
temperatures can be determined for all other points where the film 
densities are knoivn. 


Figure 5 shows the surface temperature patterns, obtained from 
thermal imagery, for the Kimberly-Clark paper mill discharge into 
the Fox River at 1:00 PM on September 16, 1971. Similar results 
have been obtained at other sites and are being used to extend the 
ground measurements to cover the whole mixing zone. 


VI . MIXING PATTERNS 


■/q aerial observations at each of the sites pro- 

measurements for defining effluent concentration 

patterns. Measurements of velocity, tem- 
perature and dye concentration and the taking of v/ater samnles at 
various locations throughout the mixing zone, along vrith plant and 

conditions, constituted the ground observations, 
eiocities were measured with Gurley, cup- type current meters having 




Figure 5 , Surface temperature contours from thermal imagery at the 

Kimberly-Clark Paper Mili mixing zone. September 16, 1971. 
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both pulse and direct reading outputs. These meters were accurate 
to + 0.1 fps. Temperatures were measured with a Whitney underwater 
thermometer having an accuracy of + O.l^C. Rhodamine WT dye, 
utilized at some sites as a tracer of effluent mixing patterns, 
was added to the effluent discharge and sensed at various points 
with a flov;- through type Turner fluorometer. Water samples of 
200 ml for suspended solids determination were taken with a specially 
built, multiple sampling device. Suspended solids were obtained by 
standard filtering and gravimetric methods. Each of these sensing 
devices was mounted in a 14 ft Jon boat used in the ground sampling 
program and was moved vertically over the river depth at each 
sampling station. These stations were generally located on sections 
perpendicular to the axis of the effluent discharge. The boat lcca= 
tion was determined by the intersection of two, emgular transit 
sightings from known positions along the shoreline. 

To delineate concentration distributions and mixing patterns, 
the field data are plotted on horizontal and vertical cross sections 
and velocity and concentration contours are drawn. Analysis of 
these plots enables the determination of the rates of effluent spread 
and dilution, namely the variation of effluent width, flow rate, 
maximum concentration and trajectory and of turbulent diffusion 
coefficients with distance from the outfall. In the paragraphs below, 
two of these mixing characteristics (variation of effluent width and 
flow rate with distance from the outfall) are presented and dis- 
cussed for several effluent, river and climatological conditions at 
the three field sites described earlier, 

Weston Power Plant . Table I stimmarizes the effluent and river char- 
ac ter is tics for five surveys. Figures 6 and 7 show the surface 
temperature distributions on two of the surveys. 

The effluent temperature and velocity distributions reveal the 
two stages of mixing described previously in Section II. Near the 
plant (principally within the bay into which the outfall discharges) 
jet-induced mixing dominates the effluent dilution. Beyond 200 to 
300 ft from the outfall, the effects of channel shape, turbulent 
diffusion and density differences (between effluent and river) 
control the spreading and dilution of the effluent. Figures 8 and 9 
show the variation of effluent surface width, b, and flov/ rate, Q, 
with distance from the outfall for the two surveys in figures 6 and 
7, These tv/o characteristics clearly show the effects of changes in 
river, effluent and climatic conditions on effluent spreading and 
dilution, namely that decreases in Oq/Q^ and increases in wind 
reduce the rate of lateral mixing. These effects can also be seen 
in figures 6 and 7. The surface width of the effluent, as determined 
by the color change between the effluent and river from 35 mm color 
aerial photographs, is in good agreement with the boat measurements, 
except for the July 15, 1970, survey when strong winds, and thus 
increased vertical mixing, made identification of the boundary from 
color change difficult. 



Table I. Characteristics at Weston Power Plant for Field Surveys 


Date 

Q^f 

cfs 

®r' 

cfs 

d, 

ft 

W, 

ft 

°C 

Qo/«r 

V /V 
o' r 

^2' 

ft X 10~^ 

Climate 

7/26/69 

187 

1800 

4.87 

533 

10.6 

0,10 

3.7 

17.4 

10.7 mph wind 

25® air, 
cloudy 

8/2/69 

187 

1600 

4.01 

518 

9,9 

0.12 

4.7 

14,1 

8 . 3 mph wind , 

25®C air, 
sunny 

7/15/70 

187 

2750 

3.71 

514 

11.3 

0.07 

2.4 

54,0 

14,0 mph wind 

22°C air, 
cloudy 

8/26/70 

187 

1800 

4.07 

503 

11.9 

o 

o 

3,6 

19,6 

5 , 2 mph wind , 

28‘^C air, 
sunny 

7/22/71 

187 

2950 

5,04 

465 

12.9 

0.06 

4.1 

47.1 

12,4 mph wind 


25^C air, 
cloudy 
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Figure 6. Surface temperature contou>"S for the Weston Power Plant on July 15, 1970 
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Figure 7. Surface temperature contours for the Weston Power Plant on August 26, 1970. 
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Fig, 8. Normalized distance from outfall versus mixing zone width and 
discharge for July 15, 1970, at Weston Power Plant. 
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Fig. 9. Normalized distance from outfall versus mixing zone width and 
discharge for August 26 » 1970, at the Weston Power Plant. 
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The vertical temperature distributions showed intense verti- 
cal mixing (vertical isotherms) on some surveys and pronounced 
tilting of isotherms on others. The vertical distributions are 
very sensitive to the intensity of the wind and secondarily to 
Qo/Qr weak winds of less than 5-7 mph, vertical 

turbulent mixing v;as suppressed by density differences due to 
temperature, which resulted in increased lateral spreading and 
inclined isotherms. For strong winds of greater than 10-12 mph, 
vertical turbulent mixing generated vertical isotherms and a 
slower rate of lateral spreading. These effects can also be seen 
on figures 6 and 7. 

Estimated lengths of the Tnjxing zone for each survey# obtained 
from a two-dimensional, steady-state model of the temperature 
distribution (including surface heat loss and channel boundaries) , 
are shown in Table I. These results show that the mixing zone 
length varies from 3 to 11 miles along the river and depends upon 
Qo/Qr wind speed , in addition to the river velocity and depth 
and lateral diffusion coefficient. In accordance with the defini- 
tion of the mixing zone, presented earlier, the temperature on the 
mixing zone boundary is given by ATp(Q^/Q ) . Extrapolation of the 
measured maximum surface temperature to this value indicates that 
the predicted mixing zone lengths in Table I are one order of 
magnitude too short. This significant discrepancy may result from 
the limited extent of the measurements, the three-dimensional 
nature of the actual temperature distribution, bends or other 
changes in channel geometry, and/or to errors in estimating the 
model parameters (diffusion coefficients, channel geometry and 
surface heat losses) . 

Waukesha STP . Table II summarizes river and effluent characteristics 
for five surveys. Figures 10 and 11 show the effluent concentration 
distribution for two surveys, as determined from temperature 
measurements. 

As the river is shallow, the effluent mixes uniformly over 
the depth due to turbulence generated by bottom frictional resis- 
tance and wind shear. Meteorological conditions, in general, appear 
to have little effect on currents or spreading patterns. Up to 
30 to 50 ft downstream of the outfall, the effluent spreads rapidly 
across the river (due to the initial momentum of the discharge 
relative to the river and the angle of discharge) to a width, b^, 
which increases linearly with Qq/Qj.. Beyond this initial region, 
and over the region of measurements, the effluent spread is con- 
trolled by lateral turbulent mixing (and to some extent, channel 
geometry) such that 

b + /2D^ x/V (6) 

where b is effluent width, Dy is lateral turbulent diffusion coef- 
ficient, X is distance from outfall and is river velocity. As 



Table II. Characteristics at Waukesha STP for Field Surveys 


Date 

°o' 

cfs 

cfs 

d, 

ft 

W, 

ft 

Q /Q 

b^/w 

V /V 
o' r 

] 

ft 

Climate 

8/9/69 

14.7 

62.1 

0.77 

87 

0.24 

0.30 

3.2 

1.6 

X 10® 

10 mph wind, 

20°C air, 
drizzle 

8/16/69 

14.1 

60,1 

0,69 

88 

0,23 

0.38 

2.9 

1,4 

X 10® 

10 mph wind, 

24°C air, 
cloudy 

8/20/69 

17,6 

55.0 

0.71 

77 

0.32 

0.43 

3.5 

1.1 

X 10® 

12 mph wind, 

20°C air, 
sunny 

8/5/70 

15.5 

33 

0.68 

77 

0.47 

0.50 

5o0 

1.1 

X 10® 

11 mph wind, 

23°C air, 
sunny 

9/10/71 

16,3 

28 

0.66 

68 

0.58 

0.75 

5.3 

8.3 

X 10^ 

15 mph wind. 


26°C air, 
cloudy 






24 


the physical measurernents covered only 1000 to 1500 ft aownstream 
from the outfall, it is not known if equation (6) is valid fartVioi- 
dovmstream (the river channel as noted earlier is very sinuous) . 
However, 35 mm color photography indicates that up to 1 or 2 miles 
dovmstream the effluent discharge remains concentrated on one bank 
of the river, and that the effluent boundary determined photograph-' 
ically corresponds to half of the effluent width by ground measure- 
ments. The flow rate within the effluent boundary increases in 
proportion to the width, b. Initial dilution, due to entrainment, 
is very rapid up to 30 or 50 ft downstream of the outfall. 

Estimates of the length of the mixing zone, obtained from a 
two-dimensional, steady-state model of the concentration distribu- 
tion (including channel boundaries) , are shown in Table II. Based 
upon these results, the mixing zone extends from Id to 30 miles along 
the river and is primarily dependent upon Qq/Q^, that is, decreasing 
as Qo/Or increases. Lateral mixing due to turbulent diffusion is 
very slov7, and a lateral diffusion coefficient of 0.23 u*d* is in 
reasonable agreement with the measured distributions. Because the 
estimated lengths of the mixing zone are very long and the measure- 
ments only covered about one percent of this distance, the validity 
of these estimates needs to be evaluated? however, aerial photographs, 
as noted earlier, show that lateral mixing is very slow. 

Kimberly-Clark Paper Mill. Table III summarizes river and effluent 
character i'^ics for five surveys. Figure 12 shows the effluent 
surface temperature distribution for one of these surveys. Figure 
5, presented earlier, shows the surface temperature distribution, 
obtained from aerial thermal imagery, on another survey. 

The effluent mixing patterns, shown in figures 5 and 12 are 
complex and result from various physical mechanisms. Near the 
plant, the 21 outfall discharges blend together and mix by entrain- 
ment of river water. Beyond 200 to 300 ft from the outfall turbu- 
lent diffusion and channel geometry and variations thereof, namely 
channel shape, alignment and depth, govern the mixing. Figure 13 
shows the variation of effluent surface width, b, and flow rate, Q, 
with distance from the outfall for the survey in figure 12 and one 
other survey. Note that these mixing characteristics are very sim- 
ilar for both surveys as the effluent and river conditions are 
nearly the same as shown in Table III, Up to about 1500 or 2000 ft 
from the outfall, b and Q increase approximately as except 

for minor variations associated with the curving bank boundary, and 
beyond this point there is rapid lateral spreading due to a radical 
change in channel alignment, that is, a sharp bend and the increased’ 
depths due to the downstream dam. The temperature and suspended 
solids distributions downstream of this bend show the development 
of large scale horizontal eddy motion as well as secondary flow. 
Further, although the temperature and suspended solids distributions 
show similar trends, their distributions, as discussed below, are 
not similar. Although b x^/^ for x < 1500 to 2000 ft, the rate of 



Table II I o Characteristics at Kimberly*-Clark Paper Mill for Field Surveys 


Date 

«o' 
cf s 

Qr' 

cf s 

d, 

ft 

w, 

ft 

^r' 

°c 

°C 

Total Solids 
Discharged, ** 

Ib/day X lO"^ 

Dy/'l*<3* 

. Local 
Climate 

8/11/71 

18.1 

1824 

8 

700 

24.0 

15.9 

40.4 

— 

10 mph wind, 

18^C air, 
cloudy 

8/27/71 

17.9 

1830 

8 

700 

20.4 

15.0 

44.4 


8 . 6 mph wind 

16®C air, 
sunny 

9/16/71 

16.9 

1640 

8 

700 

21.2 

14,8 

42.2 


7 . 3 mph wind 

12^C air, 
cloudy 

7/18/72 

16.9 

1781 

8 

700 

24.2 

00 

. 

l-l 

45.5 


6 mph wind, 

21°C air, 
light clouds 

7/19/72 

15.5 

1822 

8 

700 

24.0 

14.0 

62.5 

0.18 

12 mph vrind, 


22*^C air. 


sunny 


** 80-85 percent mineral and 15-20 percent fibre. 
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Figure 12. Relative effluent concentration contours for the 
Kimbefly-Clark Paper Mill on August 27, 1971. 
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surface spread of the effluent is 2 to 4 times that due to lateral 
turbulent diffusion. The photographic surface boundary of the 
effluent for the August 27, 1971, survey (see figure 12) is in 
good agreement with ground measurements beyond about 1000 ft from 
the outfall. The discrepancy with the ground measurements close 
to the outfall may reflect changes in plant operations or meander 
ing (large scale eddies) of the effluent boundary which is 
averaged out in the ground observations. 

The vertical distributions of temperature and suspended solids 
up to 1000 ft from the outfall show that vertical mixing is not 
sufficient to produce uniform concentrations over the depth- The 
isotherms are inclined away from shore and there are subsurface 
pockets of high suspended solids concentrations. Beyond thxs point, 
vertical mixing is intense enough to produce vertical isotherms an 
the suspended solids show slightly greater concentrations near the 
bottom. As noted earlier, the distributions of temperature and 
suspended solids are not similar. The differences in these distri- 
butions are probably due to the v/eight and resultant settling of 
the solids, forming a sludge blanket on the bottom, and to errors 
of 10 ppm in the sampling of and analysis for suspended solids. 
Finally, the only effect noted of the meteorological conditions was 
an increase in the vertical mixing in response to strong winds 
greater than 10 to 12 mph. 

No direct observations of mixing zone lengths vrere made as the 
measurements only extended to the Little Chute Dam approximately 
0.75 miles downstream of the outfall. Further, estimates of mixing 
zone lengths based upon simple mathematical models are inadequate, 
because of the effects of the bend and dam. Extrapolation of the 
Q versus x line on figure 13 indicates that 2 to 3 miles of river 
are required to dilute the effluent to the average concentration 
of the mixing zone boundary. Suspended solids concentrations of 
about 10 ppm and temperatures of about 24.5®C at the Little Chute 
Dam on the July 18, 1972 survey indicate that additional dilution 
is required. For this survey the mixing zone boundary would be 
given by suspended solids concentrations of 4 to 5 ppm and temper- 
atures of 24,3°Cr 
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NOMENCLATURE 


A,B,C, = Parameters 

* Film type, handling and processing parameter (n » 0,1, 2, 3) 
b « Width of mixing zone in surficial layer 

b » Width of mixing zone in surficial layer at the end of stage 

one (in Eq. 6, P, 19; b^ « outfaVl width otherwise) 

Hq * outfall depth 

d * Average river depth 

d* = Local river depth 

D « Image optical density 

D ® Depth averaged lateral turbulent diffusion coefficient 
E « Film exposure 

« Scanner output proportional to emitted radiation 
Kn » Parameter (n = 1,2, 3, 4) 

L.„, = Estimated length of mixing zone from a 2-D Mathematic Model 

Q = Total volumetric discharge rate of the mixing zone 

= Effluent volumetric discharge rate at outfall 
« River volumetric discharge rate just downstream from outfall 
R = Water sample reflectance 

S e Concentration of suspended solids 

T « Temperature of the surficial layer of mixing zone 

aTq ® Effluent temperature above that of the River at the outfall 

T^ » Temperature of the ambient river 

u^^ » Friction velocity or = 0,2 
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« Average effluent efflux velocity Into river at outfall 
« Average river velocity for d , and W, 

W = Average river width 

X = Horizontal distance from outfall 

Z = Relative log E 

a>$ = Model parameters 

p * Density of ambient river 
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